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NOMENCLATURE 


a = wave  amplitude 

Ac  = cushion  planform  area 

(A.).  = area  of  jupe  (cell)  bottom  opening 

J * 

A^  = leakage  area 

I(AL)i  = area  of  orifice  between  loop  and  vertical  cell 

Ag  = equilibrium  leakage  area 

b = beam 

C = orifice  coefficient 

n 

G = cushion  leakage  fraction  of  gap  height 

h^  = depth  below  origin  of  fully  extended  skirt  at  ith  jupe 

I = craft  moment  of  inertia  about  the  x-axix 

I = craft  moment  of  inertia  about  the  y-axis 

K = roll  moment 

k = wavenumber 

k^,k2  = increase  in  cell  width  on  inside,  outside  per  unit  change  in  height 

i = cushion  length 


M = pitch  moment 

m = mass  of  craft 

m = mass  of  air  in  cushion 

c 

P = cushion  gauge  pressure 

PL  = loop  pressure 

Q = volume  rate  of  flow 

Q0  1 2 = fan  performance  parameters 

AS.j  = length  of  skirt  seal  around  cushion  periphery 

T = effective  width  of  skirt  jupes 
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Uq  = craft  speed 

U = unit  step  function 

Vc  = cushion  volume 

W = craft  weight 

xc  = x coordinate  of  cushion  centroid 

zc  = depth  of  top  of  cushion  from  origin 

Zg  = heave  coordinate 

Z = force  on  craft  in  the  z-direction 

i.  L. 

B.j  = angle  of  outward  normal  at  the  itn  cell  with  respect  to  the  x-axis 

y = angle  of  wave  propagation  direction  with  respect  to  the  x-axis 

yc  = ratio  of  specific  heats  (=  Cp/cy) 

n = wave  height 

0 = pitch  angle 

x = wavelength 

p = density 

<P  = roll  angle 

= encounter  frequency 
w = wave  frequency 

Subscripts 

a refers  to  atmosphere 

1 refers  to  ith  seal  (jupe) 

j refers  to  jth  subcushion 

L refers  to  loop  or  leakage  flow 

0 refers  to  equilibrium  conditions 
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I.  Introduction 


The  purpose  of  this  report  is  to  present  the  results  of  GHR  contract  N00014- 
77-C-0063  in  which  mathematical  models  and  computer  simulations  were  developed 
to  predict  the  motion  of  air  cushion  vehicles  (ACV's)  in  random  seas.  The  craft 
are  considered  to  be  traveling  at  constant  speed  through  a sp2cified  random  sea 
of  arbitrary  heading.  The  computer  simulations  predict  the  resulting  pitch, 
heave  and  roll  response  of  the  craft.  Numerical  models  are  developed  for  two 
types  of  vehicles.  In  the  first  version,  the  craft  has  a single  cushion,  peripheral 
cell-stabilized  support  system;  in  the  second,  a divided  cushion  with  a skirt 
consisting  of  finger  seals  is  used.  The  computer  programs  are  applied  to  two 
specific  craft-the  JEFF(A)  and  the  JEFF(B)-as  particular  examples,  respectively, 
of  each  type. 

Previous  work  on  the  theoretical  prediction  of  the  seakeeping  of  air  cushion 
vehicles  which  has  been  reported  in  the  literature  includes  an  analysis  of  the 
vertical  plane  dynamics  of  SES  craft  by  Kaplan  and  Davis  (1).  They  presented  a 
simplified  approach  to  the  problem  enabling  some  discussion  of  natural  frequencies 
and  damping  ratios  to  be  made  before  resorting  to  computer  simulation.  Lavis, 
Bartholomew  and  Jones  (2)  considered  the  response  of  a peripheral  cell  air 
cushion  vehicle  in  random  seas  using  a model  which  required  empirical  data  for 
the  roll  and  pitch  stiffness  and  damping.  In  addition,  they  discussed  atmospheric 
scaling  and  the  consequent  problem  of  lift  system  parameter  distortion  in  scale 
modesl  by  using  an  analysis  of  the  flow  equations  associated  with  heave  motion. 
Linear  models  for  hovercraft  motion  in  regular  waves  were  also  developed  by 
Reynolds  (3)  and  Reynolds,  West  and  Brooks  (4)  while  a nonlinear  analytical 
model  for  the  heave  motion  of  ACV's  in  waves  has  been  proposed  by  Lebel  and 
Swift  (5).  A nonlinear  numerical  model  for  the  coupled  pitch  and  heave  motion 
of  an  air  cushion  vehicle  in  regular  waves  was  developed  by  Doctors  (6).  The 
craft  he  considered  was  of  the  divided  cushion  type  in  which  the  main  cushion  is 
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divided  into  subcushions  by  flexible  longitudinal  and  transverse  stability 
keels.  More  recently.  Doctors  has  extended  his  analysis  of  divided  cushion 
craft  to  include  the  effects  of  hydrodynamic  influence  and  compressibility 
(7,8).  Carrier,  Lundblad  and  Swift  (9)  on  the  other  hand,  have  presented  the 
results  of  a non-linear  analysis  for  the  coupled  pitch  and  heave  motion  for  both 
a divided  cushion  craft  and  a venicle  of  the  single  cushion,  peripheral  cell 
type.  More  recently,  Carrier,  Magnuson  and  Swift  (10)  have  compared  the  predictions 
of  computer  simulations  with  seakeeping  experiments  made  in  a towing  tank.  The 
nonlinear  models  used  in  references  (6-10)  required  a computer-evaluated  numerical 
solution  of  the  dynamic  equations. 

In  related  experimental  investigations,  Magnuson  (11)  and  Fein,  Magnuson 
and  Moran  (12)  reported  full  scale  seakeeping  trial  data  for  the  performance  of 
a 50-ton  British  Hovercraft.  Magnuson  and  Mesalle  (13)  have  also  used  the 
results  of  model  experiments  to  develop  linear  equations  of  motion  for  the 
JEFF(B)  in  head  waves,  while  Moran,  Fein  and  Ricci  (14)  have  investigated  the 
seakeeping  of  high  length- to- beam  ratio  SES  craft. 

In  this  report  a non-linear  model  for  pitch,  heave  and  roll  motions  in 
waves  is  developed  from  first  principles  for  a two  types  of  ACV.  The  resulting 
dynamic  equations  are  then  solved  on  a digital  computer  using  numerical  techniques 
for  two  specific  craft. 

These  craft  have  skirt  configurations  which  represent  two  different  concepts 
in  ACV  flexible  skirt  design.  Though  in  a generalized  sense,  both  concepts  make 
use  of  the  idea  of  subdividing  the  cushion  planform  area,  the  divided  cushion 
designs  employ  a small  number  of  comparable  sized  subcushions,  while  the  single 
cushion  configuration  utilizes  a large  main  cushion  with  much  smaller  cells 
contained  within  the  surrounding  skirt. 

More  specifically,  the  divided  cushion  craft  use  transverse  and  longitudinal 
stability  keels  to  compartmentalize  the  cushion  planform  into  four  sub-cushions. 
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The  purpose  of  dividing  the  main  cushion  is  to  provide  a restoring  moment  to 
perturbations  in  pitch  and  roll.  The  cushions  are  contained  by  a flexible  skirt 
around  the  craft  periphery  consisting  of  finger  seals  which  are  open  on  the 
inside  and  are  maintained  in  position  by  the  cushion  pressure.  Each  cushion  is 
fed  individually  by  fans;  the  air  then  leaks  either  to  the  atmosphere  or  beneath 
a stability  keel  to  another  compartment.  Usually  in  craft  of  this  type,  the  air 
supply  for  each  cushion  is  fed  through  a peripheral  bag  which  functions  as  an 
intermediate  plenum.  The  effects  of  the  plenum  may  be  taken  into  account  in  the 
model  by  adjusting  the  pressure  - flow  relation  for  the  cushion  feed.  The 
numerical  model  will  be  run  at  small  wave  amplitudes  so  that  no  wave  impact  on 
the  peripheral  bag  occurs. 

The  single  cushion  craft,  on  the  other  hand,  is  essentially  an  extension  of 
the  peripheral  jet  concept  by  use  of  flexible  trunks  to  maintain  adequate  clearance. 
In  this  configuration,  a single  cushion  is  contained  by  a flexible  skirt  made  up 
of  a series  of  cells.  Air  is  fed  through  orifices  to  the  cells  or  jupes  from  a 


loop  plenum  located  above  the  jupes  and  extending  around  the  craft  periphery. 

The  air  supply  inflates  the  cell  and  then  escapes  through  the  open  bottom  forming 
a vertical  jet.  The  loop  itself,  as  well  as  the  main  cushion,  are  each  supplied 
from  a separate  fan  system.  The  cellular  skirt  system  provides  the  same  static 
stability  function  as  do  the  stability  keels  of  the  divided  cushion  craft.  In 
addition,  the  downward  jet  flow  beneath  the  skirt  around  the  craft  periphery 
forms  a "curtain"  effect  minimizing  leakage  from  the  main  cushion. 

In  this  report,  a mathematical  model  for  each  type  of  craft  is  developed 
separately  in  the  next  two  sections.  This  is  followed  by  a section  describing 
the  synthesis  of  the  sea  state  excitation.  In  the  last  part,  the  computer 
programs  obtained  by  solving  the  dynamic  equations  numerically  are  applied  to 
the  JEFF(A)  and  the  JEFF(B).  The  computer  program  documentation,  however,  is 
contained  in  a separate  report. 


II.  Peripheral  Cell -Stabilized  ACV 

The  position  and  orientation  of  the  ACV  is  given  by  the  heave  coordinate  Zq 
(positive  downward),  pitch  angle  e,  and  roll  angle  <(>.  An  x,  y,  z coordinate 
system,  fixed  with  respect  to  the  craft  is  located  so  that  its  origin  is  at  the 
center  of  mass  (see  Figure  1).  The  craft,  traveling  at  constant  speed  Uq, 
encounters  random  seas  in  which  the  surface  height,  with  respect  to  the  mean  ele- 
vation, is  given  by 


n = n(x,y,t) 


The  equations  of  motion  governing  the  craft  response  are: 


m z0  ^cushion+^seals+W 


(1) 


*y  9 ^cushion+^seals 
*x  * ^cushion+^seals 


(2) 


The  cushion  and  seal  forces  and  moments  will  be  specified  in  terms  of  craft  position 
and  cushion  and  cell  pressures. 

The  main  cushion  pressure  must  be  such  that  the  rate  of  change  of  the  cushion 
mass  equals  the  net  rate  of  mass  flow  into  the  cushion.  This  is  expressed  by: 


dm 

c 

IT 


pin^in  " °out^out 


(3) 


It  will  be  assumed  that  the  air  density  while  entering  and  leaving  the  cushion  is 
approximately  the  atmospheric  air  density. 


in  pout 


= p. 


(4) 


The  flow  into  the  cushion  is  due  to  the  fan  air  supply  and  is  expressed  in  terms 
of  the  cushion  pressure 


f 


(5) 


^in  = Q0  + Q1P  + Q2P 


where  the  coefficients  are  determined  empirically.  Equation  (5)  approximates  the 
"fan  curve"  for  the  supply  fans. 

The  flow  out  of  the  main  cushion  is  due  to  leakage  between  the  skirt  and  the 
water  surface  and  is  specified  using  an  orifice  coefficient  modification  to  the 
Bernoulli  equation 


’out  * CA  \J^ 


(6) 


The  total  leakage  area  A^  consists  of  an  equilibrium  leakage  area  Ag  and  leakage 
area  due  to  the  physical  gap  between  each  element  of  the  skirt  and  the  surface, 

«Wo 


AL  = Ag  + G I H.aS. 


A0 


craft 

periphery 


Hi  = ^VVV'V'^H("VVV"V_r,i)’  ni  = 


(7) 


The  static  area  AQ  accounts  for  leakage  area  due  to  the  cushion  pressure  in- 
duced free  surface  deformation  at  equilibrium. 

The  gap  height  at  the  ith  cell,  Hi  is,  however,  not  all  available  for  main 
cushion  leakage.  Air  from  the  underfed  peripheral  jets  directed  downward  beneath 
the  cells  must  also  escape  under  the  outer  hem  of  the  skirt.  That  is,  the  actual 
flow  through  the  physical  gap  between  the  outer  skirt  and  the  water  surface  consists 
of  flow  from  the  main  cushion  and  deflected  jet  flow  from  the  cellular  skirt.  The 
effective  gap  height  for  the  main  cushion  is  then  a fraction  G of  the  physical 
gap  height.  The  proportion  of  the  gap  alloted  to  main  cushion  flow  will  be 
approximated  by  the  fraction  known  to  exist  at  equilibrium.  This  would  be  equal 
to  the  ratio  of  equilibrium  cushion  supply  to  the  total  fan  feed  at  equilibrium. 

The  rate  of  change  of  the  mass  of  air  in  the  main  cushion  is  expanded  by 
writing 
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drn  d\/  . 

c _ „ c u dp 

W ’ p dT  + Vc  3t 

The  cushion  volume  depends  on  the  craft  position  relative  to  the  waves: 


Vc  = Ac(~zO+xc0-zc)  - ||  n(x’*’t)dA 


C r.hion 
Planform 

The  first  term  in  the  above  expression  is  due  to  changes  in  average  height  of 
the  cushion,  while  the  second  term  represents  the  wave  "pumping"  effect. 

Compressibility  effects  of  the  cushion  air  mass  are  taken  into  account  by  the 
last  term  in  Equation  (8).  The  cushion  density  is  related  to  cushion  pressure  by 
assuming  that  the  thermodynamic  processes  are  reversible  and  adiabatic.  The  rate 
of  density  change  may  then  be  expressed  by 


dp  _ 1 1-y  /Pa  x dP 

art  — 7 p at 


Thus,  Equations  (3)  - (10),  which  represent  conservation  of  mass  for  the  main 
cushion  air  flow,  relate  the  cushion  pressure  to  craft  position,  surface  wave  motion 
and  the  rates  of  motion  of  the  vehicle. 

The  cushion  forces  on  the  craft  are  the  resultants  of  the  cushion  air  pressure 
acting  over  the  area  of  the  cushion  supported  externally.  The  z-direction  force  is 


^cushion  -PAC’ 


while  the  pitch  and  roll  moments  are 


'"’cushion  ’ l '-cos8i’ptr'i-xi6tJ'i*’lSi(-z0-'’i>-xc(Z>cusMl 


'^’cushion  ' l (^"S^PCnj-x^y^aSif-Zg-n,) 


respectively. 
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The  seal  forces  will  be  specified  by  considering  each  jupe  or  cell  in  a manner 
roughly  analogous  to  the  treatment  of  the  main  cushion.  The  open  cells  themselves 
may  be  thought  of  as  a series  of  mini-cushions  around  the  bow  and  sides  of  the 
craft.  For  the  i*  cushion,  conservation  of  mass  is  required  for  the  air 
flow  through  the  jupe.  This  would  be  a statement  that  the  volume  rate  of  flow 
into  the  cell  from  the  loop  plenum  must  equal  the  flow  out  the  cell  bottom  opening 
plus  the  rate  of  cell  volume  increase.  In  equation  form,  this  becomes 

wv- 

+ T AS. 

( - z 0+x i © ~y i - n i )uSp(zo"xie+yi<f>+hi+rii ) O4) 

i.  L_ 

where  (AQ)..  is  the  equilibrium  leakage  area  under  the  iin  jupe.  Since  the  equilibrium 
jupe  pressure  is  the  same  as  for  the  main  cushion,  (AQ)^  can  be  expressed  in  verms 
of  the  cell's  equilibrium  flow  rate  by 


For  the  closed  jupes  across  the  stern,  the  second  term  in  Equation  (14)  is  omitted 
since  there  is  no  escape  flow  out  the  bottom. 

Though  the  loop  plenum  is  fed  separately  from  the  main  cushion,  the  gap  exit 
flow  areas  for  the  cushion  and  the  loop-skirt  system  are  proportional.  The  flow 
out  from  the  two  systems  is  also  nearly  proportional  since  the  pressure  of  the  cell 
is  usually  very  close  to  the  cushion  pressure.  Because  the  total  flow  through  the 
loop-skirt  system  is  nearly  proportional  to  flow  through  the  cushion,  the  loop 
pressure,  PL  a P.  The  constant  of  proportionality  will  be  taken  as  the  ratio  of 
the  respective  pressures  at  equilibrium.  With  this  approximation  and  conservation 


Cn<\>i 


^77 


sgn  [PL-(Pj)f]  - cn[(A0)+(i-G)as1(- 


-hi-n1  )]us  [(A0)1+O-G)aSi(-z0+x1e-y1+-h1-ri1)] 


2<Vi 


of  mass,  the  pressure  within  each  cell  may  be  related  to  cushion  pressure,  craft 
motion  and  wave  motion. 


The  additional  vertical  force  (above  that  already  included  in  the  main 
cushion  forcing  term)  on  the  itfl  seal,  Z.,  is  expressed  in  terms  of  the  cell 
pressure  by 


(16) 

(17) 

III.  Divided  Cushion  AC V 

The  motion  of  a divided  cushion  ACV  in  the  arbitrary  wave  field  (specified 
formally  by  Equation  (1))  is  governed  by  the  same  equations  of  motion.  Equations 
(2)  as  the  single  cushion  craft.  The  cushion  and  seal  forcing  terms  as  well  as 
sub-cushion  conservation  of  mass  expressions  will,  of  course,  be  different. 

Let  subscript  j represent  either  I,  II,  III  or  IV  corresponding  to  the  sub- 
cushions shown  in  Figure  2.  Then  the  conservation  of  mass  expression  for  the 
jth  cushion  is  given  by 


= - [(Pj)i-P](Aj)i+{-[(Pj)i-P](k1).[z0-xie+y.<j.+h.+ni]- 
(Pj^-(k2^,[zO-xi0+yi*+hi+ni]}ASiUsp(zO'xi0+yi4,+hi+ni) 


The  resultant  seal  forces  and  moments  are  then: 


'Z’seal  * l V '"’seal  ' J'W  <K’seal  * I'*i’zl 


at  <mc>  = pa[(Qm).  ‘ 'Qout’ 

■ S'  K> 

■ pj  at  «'c>  + <»e>  at  -j  <’8> 

J J 

In  the  craft  considered,  each  subcushion  is  supplied  air  separately  by  fans. 

i.  L. 

The  air  supply  for  the  j cushion  is  expanded  to  second  order  in  terms  of  the 
corresponding  cushion  pressure: 
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= Coj  + Clj  Pj  + C2j  Pj‘ 


where  the  empirical  constants  2 3j  include  the  effects  of  ducting  from  fan  to 


cushion. 


The  flow  out  of  the  j cushion  due  to  leakage  between  skirt  and  water  surface 


Is  expressed  by 


5°ut’j  ' C"  (Al,j  * c"  Jl 


V.  sgr<pj-pm) 


The  leakage  area  (A,  ) represents  gap  openings  between  the  jtn  cushion  and  the 

j 

atmosphere. 


(A,  ) • (A.)  + ( l H.  AS.) 
j j 1 j 

External 
Periphery 


where  the  summation  is  around  that  portion  of  the  skirt  separating  cushion  j from 


the  atmosphere.  The  area  AmJ.  is  the  leakage  area  between  cushion  m and  cushion 


j and  is  expressed  by 


Vi  ■ <£  H1  iSi>, 


in  which  the  summation  is  across  the  skirt  separating  cushion  m from  cushion  j. 


The  effective  static  leakage  area  for  the  cushion,  (AQ)j»  is  given  by 


Vat  V 


and  accounts  for  free  surface  compliance  effects  on  the  leakage  area  while  the  craft  Is  In 

w 

equilibrium.  The  pressure  term  Rj  is  the  static  gauge  pressure  of  cushion 

j.  (p0)j»  and 
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when  j = 1, 1 1 


<_Xc,III 


Rj  ■ 


(24) 


(xc), 


2L(-Xc),ii  ♦ IW 


when  j = III, IV 


To  simplify  the  analysis,  the  cross  flow  terms  are  linearized  using  an  equivalent 
linearization  procedure  with  the  result  that 

Wi  - U 


'n  ^lj  f 


ip  _p  ) = nC,n  V)  (P  _P  ) 

sgn  <Pj  V XT?  1 j 


(25) 


mj  "a 


in  which  V . is  a characteristic  crossflow  velocity  and  is  taken  to  be  \[ 
mj  1 

1L 

The  reversible,  adiabatic  assumption  applied  to  the  j cushion  yields 

d .1  ^/"aV'd  P,, 

l C jj 

while  the  volume  of  cushion  j,  (V  ) , is  given  by 

C j 

(VJ  = A .[-z  + (xj  e-  (y_)  ♦-  zr] 


(26) 


CJ  ^ 0 


C j 


Plan 


(27) 


n(x,y,t)dA 
orm 


Cushion  j 

The  cushion  forces  and  moments  are  the  resultants  of  the  cushion  pressure  acting 
over  externally  supported  areas: 

IV 


^cushion  " ^ 

IV 

(M)cushton  I,  [ { (-“s  SifjOVV-V’' 


periphery 
cushion  j 


+ pj  Aj  (xc>j 
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J 


IV/- 

'"'cushion  = jljJ  ' $ <s1"  6i>  Pj  '"rx1,-»1*,4S1 


periphery 
cushion  j 


("zo‘ni"j  - Pj  flj  (*c'j] 


(28) 


The  seal  forces  are  calculated  assuming  that  the  seal  fabric  has  no  Inertia 
so  that  the  hydrodynamic  pressure  on  the  seal  equals  the  air  pressure  contained  by  the 
skirt;  thus: 

IV  % „ 

(Z)seals  ' I,  - pj  t J MVV+V*V"l'VVx19+V+’’,>4Si:lj 

external 

periphery 


IV 


M 


seals  Pj  [ ] xikl^o"x1e+y1*+h1+ni^Usp^2o"x10+yi*+h1+n1^S1^j 

external 
periphery 

IV  „ , 

seals  = Pj  [ l (-y1)ki(zo-x10+y1*+h1+n1)Usp(Vxie+y1+h1+n1)AS1]j 

external 

periphery 


IV.  Wave  Surface 

The  wave  surface  height  with  respect  to  the  mean  elevation  is  specified  in  the 
form  of  a finite  Fourier  sine  series: 


(29) 


N 


m=l 


n(x,y,t)  = l am  sin  [x  cos  ym+y  sin  yj  - (u>e)  t + 


(A* 


~e'  ’ Tm 

m 


j 


(30) 


By  using  this  general  Fourier  representation,  any  arbitrary  sea  surface  may  be 
generated  as  a superposition  of  small -amplitude,  deep  water  waves.  For  the  m^ 
harmonic  constituent,  the  encounter  frequency,  wavenumber  and  phase  velocity  are  related 
by  the  following  expressions: 
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— rm 


(u_)_  = o»_  - “m  Uo  cos  y, 
' e m m — — 1 


m 


(3,) 

Each  wavelet  is  defined  when  its  amplitude  a^,  heading  Ym»  frequency  a>m,  and 
phase  angle  are  given.  This  method  offers  considerable  flexibility  In  the 
specification  of  the  sea  surface.  Single  frequency  components  may  be  used  as  well 
as  super  positions  of  regular  waves. 

For  a fully  developed,  random  sea  the  Plerson-Moskowitz  spectrum  can  be  applied 
to  obtain  the  amplitudes.  In  this  approach,  the  sea  Is  assumed  to  consist  of  long- 
crested  waves  propagating  at  a single  heading  y.  The  random  sea  is  synthesized 
by  superimposing  N harmonic  wave  components,  as  in  Equation  (30),  in  which  the  com- 
ponent phases,  <f>m,  are  chosen  as  random  numbers  with  equal  probability  between  0 
and  2tt.  The  range  of  frequencies  for  which  there  is  appreciable  wave  energy  is 
divided  into  N bands  of  width  Aa>  and  center  frequency  a>m.  Then  the  mth  wave 
component  will  have  wave  frequency  u>m  and  amplitude  a^  given  by 

am  = ¥ S ^m,Aw  (32) 


where 


a 9 

s<“m>  ■ -v 


exp 


m 


f" 


(33) 


is  the  Pierson-Moskowitz  wave  energy  spectrum.  The  parameter  a=8. 10x10 
and  3 = .74,  while  V is  the  speed  of  the  wind  generating  the  sea. 


-3 


V.  Numerical  Model  Application 

The  dynamic  equations  obtained  from  the  theroetical  model  are  solved  numer- 
ically using  a digital  computer.  After  experimenting  with  several  algorithms 
for  solving  differential  equations,  it  was  found  that  a polynomial  extrapolation 
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method  (taken  for  Gear  (15))  was  the  best  means  for  extrapolating  time  dependent 
variables  from  one  time  step  to  the  next.  The  details  of  the  computational  pro- 
cedure, including  complete  listings  of  the  programs,  are  contained  in  the  "Computer 


Program  Manual  for  ACV  Simulations." 

In  the  simulations,  the  craft  trajectories  are  calculated  as  the  vehicle 
runs  at  constant  speed  and  course  through  a specified  wave  field.  The  time 
response  of  dynamic  variables  is  then  analyzed  to  determine  averages  and  root 
mean  squares.  Also  the  first  harmonic  magnitude  and  phase  is  calculated  at  each 
encounter  frequency  used  in  specifying  the  wave  field.  These  results  are  normalized 
by  dividing  the  magnitude  of  the  response  variable's  first  harmonic  by  the 
corresponding  quantity  in  the  wave  field  excitation. 

The  computer  programs  were  applied  to  two  specific  craft  - a single  cushion 
craft,  the  JEFF(A),  and  a divided  cushion  craft,  the  JEFF(B).  The  craft  particulars 
for  each  example,  as  used  in  the  programs,  are  given  in  Tables  1 and  2 respectively. 

Before  the  programs  were  run  using  random  seas,  they  were  tested  using  a 
single  frequency  component  wave  excitation.  The  procedure  was  to  run  the  craft 
at  constant  speed  through  regular  waves  until  a steady  state  response  was  achieved. 
This  was  repeated  for  various  speeds,  waveheights,  wave  lengths  and  wave  propagation 
directions.  The  results  were  then  presented  in  the  form  of  frequency  response 
plots.  Figures  (3)  and  (4)  show  a frequency  response  for,  respectively,  the 
JEFF(A)  and  the  JEFF(B)  obtained  while  the  craft  were  encountering  head  seas. 

The  JEFF(A),  in  particular,  was  extensively  studied  in  this  manner  since  towing 
tank  experimental  data  was  available  for  comparison.  The  results  of  this  study 
were  presented  in  reference  (10). 

The  computer  programs  were  then  tested  using  random  sea  input.  The  data 
for  two  runs,  one  for  each  craft,  are  presented  in  Tables  3 and  4.  It  was  found 
that  the  computer  programs  operated  successfully  when  the  wave  input  was  changed 
from  regular  waves  to  the  superpositions  of  regular  waves  used  to  synthesize 
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Table  1.  Physical  Parameters  for  JEFF(A) 


W 

I, 

I, 

i 

b 

h 


k 

k 


1 

2 


T 


% 

<>1 

q2 


(Vi 

(Vi 


= 320,000  lbs. 

= 4.75  x 106  ft-lb-sec2 
= 1.21  x 106  ft-lb-sec2 
= 82.2  ft. 

= 42  ft. 

= 10  ft. 

= 0 ft. 

= 5 ft. 

= +90°,  -90°  along  starboard,  port  side  of  craft; 
0 at  bow;  180°  at  stern 

= 1 along  bow  and  sides;  4 at  stern 

= .6  along  bow  and  sides;  4 at  stern 

= 4.33  ft. 

= 8.0  ft. 

= 4200  ft3/ sec 

= -11  ft5/lb-sec 

= 0 

= 2120  lb/ ft2 
= 2.2  x 10“3  slug/ft3 
= 1.3 
= .7 
= 1/4 
=2.5  ft2 

= 20  ft2  for  bow  and  sides;  0 at  stern 
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f 


Table  2.  Physical  Parameters  for  JEFF(B) 


= 325,000  lbs. 

= 10,093  slugs 
= 4. 37x1 06  ft-lb-sec2 
= l.lxl 06  ft-lb-sec2 
= 76.4  ft. 

= 38.2  ft  for  j = I - IV 


(xc}j 

(yc>j 


= 48  ft. 

= 24  ft.,  j * I - IV 

= 9.5  ft  on  peripheral  skirt,  8.5  on  stability  keels 

= > (Xch,II  = 19,1  ft,;  ^chll.IV  = ”19,1  ft‘ 

= > ^Yc^I,IV  = 12  ft';  (VlI.III  = "12  ft‘ 


{ZC^j 

^in^j 


= 4.5  ft.,  j = I-IV 
= 2010  - 3.5  P. 

= 2010  ft3/sec,  C,.  = -3.5  ft3/sec  ( J-*-),  C?. 

1 b/f tc 

= 14.7  psi  = 2.12xl03  lb/ft2 
= 2.2xlO"3  slug/ft3 
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Normal i zed 
Response 
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10  5 \/i  1 


Figure  4 JEFF(B)  head  seas  frequency  response 
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random  seas.  Due  to  limitations  on  the  amount  of  uninterrupted  computer  time 
available  however,  it  was  not  possible  to  make  random  sea  runs  which  are  long 
enough  to  yield  completely  reliable  statistics  for  the  ACV's  random  sea  motion. 
Hence,  the  short  runs  described  in  Tables  3 and  4 should  be  interpreted  as  a 
demonstration  of  the  programs'  capabilities  only. 


Table  3.  JEFF(A)  Response  to  a Random  Sea 
Record  Length  =11.1  seconds 


Wave  Input  (Head  Sea) 

Craft  Response 

Encounter  frequency 

Ampl i tude 

Phase 

Heave 

(Normalized 

Pitch 
1st  Harmoni 

1.09 

0.027 

5.83 

0.79 

0.53 

1.40 

0.245 

4.58 

0.73 

0.94 

1.75 

0.566 

5.72 

0.95 

1.13 

2.13 

0.761 

5.32 

0.79 

0.79 

2.55 

0.807 

4.03 

0.88 

0.95 

3.00 

0.766 

0.52 

0.73 

0.90 

3.49 

0.690 

4.13 

0.60 

0.57 

4.02 

0.607 

5.38 

0.41 

0.38 

4.58 

0.528 

4.02 

0.30 

0.13 

5.18 

0.458 

2.46 

0.22 

0.11 
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Table  4.  JEFF(B)  Response  to  a Random  Sea 


Fr  = 1.15 

Record  Length  = 13.5  seconds 

Encounter  frequency 

Wave  Input 

Amplitude  Phase 

Craft  Response 

Heave  Pitch 

(Normalized  1st  Harmonic) 

1.09 

0.027 

5.83 

1.37 

1.38 

1.40 

0.245 

4.58 

1.03 

1.07 

1.75 

0.566 

5.72 

0.94 

1.01 

2.13 

0.761 

5.32 

0.90 

0.90 

2.55 

0.807 

4.03 

0.86 

0.87 

3.00 

0.766 

0.52 

0.76 

0.77 

3.49 

0.690 

4.13 

0.65 

0.71 

4.02 

0.607 

5.38 

0.65 

0.72 

4.58 

0.528 

4.02 

0.47 

0.60 

5.18 

0.458 

2.46 

0.17 

0.48 
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